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Abstract- This paper proposes a robust design scheme for the 
control of surge phenomenon appearing in a centrifugal 
compressor. It is known the surge instability in compressor might 
lead to reduce the operating efficiency and even damage the 
machinery while both of unstable compressor characteristics and 
compressor load torque are hard to exacdy measure in the 
practical application. Based on the assumption of local sector 
nonlinearity for system dynamics, a robust control law is 
proposed to cover the system nonlinearity and uncertainties while 
guaranteeing the stability of system equilibrium. This is achieved 
by using the actuation of the driving torque and closed-couple 
valve. Numerical simulations are given to demonstrate the success 
of the proposed robust controller design. 

Keywords- Robust control; Centrifugal compressor; Surge 
phenomenon; Compressor characteristic; Local sector nonlinearity 

I. Introduction 

Over the past several years, the compression systems have 
been widely applied to civil and industrial applications such as 
air conditioner, automotive engines, aerospace and thermal 
plasma torch (e.g., [1], [2]). Among those applications, the 
compressors might be operated in the high pressure rise 
regarding the performance and efficiency requirements. Surge 
and rotation stall are two basic types of instability, which are 
known to appear in axial compression systems and might 
reduce the operating efficiency and even damage the 
machinery (e.g., [3]-[7]). In order to study those instabilities, 
Greitzer derived a nondimensional axial flow compressor 
dynamical model from which surge and rotation stall can be 
predicted [3], Moore and Greitzer employed the Galerkin 
method to construct a third-order model for the study of 
rotating stall [4]. Among the existing studies, several control 
strategies have been proposed for the control of unstable 
instabilities in an axial compressor by adopting the Moore- 
Greitzer model (e.g., [5]-[7]). The axial flow compression 
system described by [4] is known to be uncontrollable (e.g., [5], 
[7]). Liaw and Abed proposed a local bifurcation analysis and 
bifurcation control to stabilize the instabilities of axial flow 
compressor [5]. In the study of [5], the existence conditions of 
surge and rotation stall were found by applying local 
bifurcation analysis. To further study the prevention of the 
unstable dynamics, Liaw and Abed employed the active 
nonlinear control to transform the subcritical bifurcation into a 
supercritical one. In addition, Chaturvedi and Bhat presented 
the output feedback backstepping controller that combined a 
high gain observer to estimate the mass flow [6]. Based on the 
dynamical characteristics of practical system, the unstable 
branch of compression characteristic is difficultly measured so 
that there exists an uncertainty term in an axial flow 
compressor. Liaw and Huang proposed a robust control 



strategy to stabilize instability via sliding mode design with 
respect to the uncertainty appearing in the axisymmetric 
characteristics maps [7]. 

The study of a centrifugal compression systemhas attracted 
considerable attention in the recent years (e.g., [8]-[13]). 
Leonessa et al proposed a lumped-parameter surge model for a 
centrifugal compressor that involves the state of spool speed, 
pressure, and mass flow rate [8]. In order to prevent the surge 
dynamics, a switching type of globally stabilizing controller 
was also applied to control the unstable dynamics in a 
centrifugal compressor. In the study of [9], Gravdahl and 
Egeland developed an approximate compressor characteristics 
which considers the losses of incidence, friction, back flow, 
clearance etc. Besides, a surge control law was designed for 
variable speed operation of centrifugal compressor via PI speed 
control [9]. Bohagen and Gravdahl presented the surge control 
scheme with respect to driving torque as a control input via 
backstepping and passivity approaches [10]. In addition, the 
centrifugal compressor might exist the uncertainty term in the 
characteristic maps, which are similar to that in axial flow 
compression system. 

For the study of centrifugal compressor with system 
uncertainty, the robust control of surge behaviour has been 
proposed in [11]-[13]. Daroogheh et al employed a high gain 
robust adaptive controller to suppress the surge dynamics of 
centrifugal compression system with variable speed [11], In 
that investigation [11], non-dimensional time varying 
disturbance rejections are considered in compressor surge 
control. The control schemes in [12] were based on the 
backstepping method to achieve global stability. In order to 
compensate the system uncertainty, robust global stability via 
the Lyapunov redesign technique was employed in such a 
circumstance [12], In addition, Shehata et al provided an active 
control in a constant speed centrifugal compression system 
which is simplified to the second order model [13], That study 
employed a different type of controller based on the specific 
requirements of performance. Among those studies, robust 
control schemes for centrifugal compressor have been 
proposed to deal with uncertainty or disturbance. However, the 
robust design proposed in those studies requires the bounds of 
system uncertainty which might not be easily obtained in 
practical applications. 

In this paper, we consider different approach by using 
sector nonlinearity approach. The modelling error and/or 
measurement error of compressor's driving torque and 
characteristic maps will be treated as system uncertainties. 
According to the study of [7], Liaw and Huang utilized an 
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"additive-type" control to the mass flow of axial flow 
compressor such that the domain of attraction for stable 
equilibrium was enlarged via sliding mode control. Motivated 
by [7], one of the main goals of this paper is to enlarge the 
domain of attraction for providing the non-local stability of 
system equilibrium in centrifugal compressor dynamics. By 
such a design, the system uncertainty will then be compensated 
for guaranteeing system stability. In this paper, we try to use 
the local sector boundary to cover the nonlinear term in the 
compressor system. The local sector nonlinearity guarantees 
that nonlinear terms of the system can be bounded by a linear 
function (e.g., [14]-[15] ). 

This paper is organized as follows. In Section II, we recall 
the third-order nonlinear lumped-parameter model for 
centrifugal compression system of which the instabilities are 
described. The robust control schemes are derived in Section 
III. Section IV reveals the numerical simulations which 
demonstrate the feasibility of the proposed control laws. 
Finally, the conclusions are summarized in Section V. 

II. Dynamical Behaviour of a Centrifugal Compressor 

A centrifugal compressor is known to be consisted of inlet 
duct, outlet duct, plenum, exit duct, valve, and compressor 
motor as depicted in Fig. 1. Detailed modelling process 
through instrumentation and experimental results of centrifugal 
compressor can be referred to (e.g., [8], [9], [16]). A nonlinear 
lumped-parameter third-order equation by using notation of [12] 
is recalled as follows: 



J 



dco 



dt 
L dm c 
A dt 

dt V, 



-(rh c -rh t ), 



(1) 
(2) 

(3) 



where m denotes spool angular speed while r < and c are the 
drive torque and compressor load torque, respectively. J 



denotes the spool moment of inertia. 



is the pressure 



downstream of the compressor. p is the pressure within the 
plenum. L is the length between the compressor and duct while 

A is a reference area. ' and c denote the mass flow rate 
through the throttle and mass flow rate at the plenum entrance, 
respectively. In addition, a denotes the inlet stagnation sonic 

V 
velocity and p is the plenum volume. Note that, forthe case of 

the spool speed of centrifugal compressor being set to a 

constant value, the system (l)-(3) is found to be equivalent to 

the model derived by Greitzer (1976) [3]. 



Compressor 



Throttle 



Air flow 




Inlet duct 



Outlet duct 



It is known that (e.g., [8], [9], [16]) that the compressor 
load torque is a function of the spool speed and the mass flow 
rate. The system (l)-(3) can be rewritten in non-dimensional 
form as given below 



a> = a[r t -T(co,(f>)] 
(f> = b[y/ c (a>,(f>)-y/] 
iy = c[0-F(y,y/)], 



(4) 
(5) 
(6) 



where ^' ' and ^ c ^ '?' are defined as the compressor 
load torque function and the compressor characteristic map, 
respectively. It is observed from equations (4)-(6) that the 
nonlinear terms mainly originate from compressor load torque 

*• '^ , characteristic maps ^ c ^ '^' and throttle opening 

characteristic function ('>'/') . The dynamics of load torque 
and characteristic map in compressor are known to depend on 

the effects of mass flow and spool speed. In addition, ('>V') 
is generally a strictly increasing function with respect to both 
of the pressure rise and the throttle opening''. Moreover, in 

practical application, we have a > " , o > , and c > ® , 
which are the non-dimensional parameters defined in (e.g., [8], 



[9], [16]). Denote 



x e =(a) e ,0 e ,y/ e ) 



an equilibrium point of 



systems (4)-(6). It is known from [12] that the equilibrium 
point e of systems (4)-(6) will be asymptotically stable if 



dy/ c (co,<f>) 



<0 



dy/ c (co,(j)) 



>0 



e <* ' but being unstable if Sf * '* e .Normally, 
when a compressor operates near the surge line, a small 
variation of system parameter might cause the compressor to 
exhibit an oscillative behaviour which is found to be related to 
compressor characteristic maps [12]. Unfortunately, the 
unstable compressor characteristics are hard to measure in 
practical applications. The curve fitting scheme is usually used 
to approximate the unstable characteristic maps. Such an 
approach has been used in [8] for obtaining the characteristic 
maps of centrifugal compressor with respect to different setting 
value of spool speed. One example is shown in Fig. 2. 

The bifurcation diagrams for system (4)-(6) with respect to 

the variation of compressor torque r < adopted from ([8], [12]) 
are shown in Fig. 3, where the notations HBand PD denote the 
Andronov-Hopf bifurcation point and period-doubling 
bifurcation point, respectively. In addition, the solid line (resp. 
open circle) denotes the stable equilibrium branch (resp. stable 
limit cycle) while dot line (resp. solid circle) denotes unstable 
equilibrium branch (resp. unstable limit cycle). It is observed 
from Fig. 3 that the system states exhibit a stable limit cycle 
when the compressor operates at the equilibrium point near the 
bifurcation point HB. The timing responses from numerical 
simulations are shown in Fig. 4 to demonstrate such a 



0.1 



0.315 



Fig. 1 A Schematic Diagram of Centrifugal Compression System 



behaviour. Here, we choose t( u ' x and Y U - J±J (i.e. 
(£ e ,^) = (0.557403,0.199432,0.4008) ) ^ initial yalue 

(£ ,<Wo) = (0.5,0.345,0.22) _ The oscillative behaviour 
depicted in Fig. 4 is known to be the so-called "surge" 
instability. 
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Fig. 2 Characteristic Maps of a Centrifugal Compression System 
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Fig. 3 Bifurcation Diagrams with Respect to the Variation of Driving Torque 

Timing response Timing response 




Fig. 4 Timing Responses andPhase Portraits 

It is shown in [17] that the variation of the driving torque 

f and/or the throttle opening ' might cause a centrifugal 
compressor to exhibit surge dynamics. In addition, the domain 
of attraction for the system equilibrium points is found in [17] 

to become smaller as the compressor drive torque T ' decreases. 

Due to the effects of uncertainty in ^ '™' and ^ c ^ '^' , 
those will result in the system running into the regime of HB. 



Thus, the surge instability might be avoided by making the 
domain of attraction large enough. In this paper, a surge 
recovery control schemes is employed to drive system 
dynamics toward a stable equilibrium by using compressor 
drive torque and the actuation of closed-couple valve. Besides, 
the proposed control scheme will also enlarge the domain of 
attraction for the stable equilibrium point, which will then keep 
the system behaviour in the stable manifold disregard the 
appearance of disturbance and/or system uncertainty. 

III. Robust Control Scheme 
In this section, we will propose robust control laws to 
stabilize the compressor dynamics with system uncertainty by 
using the actuation of compressor driving torque and close- 



coupled valve. Let 



x e ={a e , e , y/ e ) 



be the equilibrium point 



of system (4)-(6) for the given driving torque <■ te and 

throttle opening ' ~~ ^ e . Here, we consider the control inputs of 
the compressor will be the compressor driving torque and the 
close-coupled valve. 



Let 



— \ 1 ) 5 J '-t / 



With ! ~ W ~ We , *2-V-Ve 

te and 2 denotes the actuation of the 
close-coupled valve. The system (4)-(6) can then be rewritten 
as follows: 



H = V~W e «! = *■( 



x 1 =a[f(r te ,x 1 ,x 2 ) + u 1 ] (7) 

x 2 =b[C ss (x 1 ,x 2 )-x 3 + u 2 ] (8) 

x 3 =c[x 2 +g(y e ,x 3 )], (9) 
where 

f{r K ,x l ,x 2 ) = T K -T{x 1 + a> e ,x 2 +<j> e ) (10) 

C ss( x l> x 2) = V / e( x l + a>e' X 2 + 0e) ~ Ve(®e'0e) C 11 ) 

g(r r , x ,) = F(re,v / e)-F(r r , x 3 +Ve) ( 12 ) 

with the following equilibrium conditions for system (4)-( 6): 

r tt =r(£ e ,£) (13) 

<t> e =F(y e , ¥e ) (14) 

y/ e =y/ c (a> e ,(j> e ) ■ (15) 

The main goal of this paper is to design control laws by 

using control inputs 1 and 2 to compensate the surge 
instability and enlarge the domain of attraction for the 
stabilized system equilibrium Detailed are given as follows. 

A. System Uncertainty for a Centrifugal Compressor 

In the practical application, the stable portion of the 
dynamical characteristics for a centrifugal compressor can 
easily be measured from the experimentation. On the contrary, 
the unstable behaviour might be difficult to obtain from the 
practical compression system Due to the fact of the inexact 
measurement, it is unavoidable to have the appearance of the 
uncertain terms in the compressor dynamics for practical usage. 
In addition to compressor characteristic maps, the uncertainty 
in compressor load torque might also emerge from the 
fluctuations of the unstable spool angular speed and mass flow 
rate. Thus, both compressor load torque and characteristic 
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maps are required to consider the existence of uncertain terms 
in the robust controller design. Let the compressor torque and 
characteristic map are given as (e.g., [12]) 

f(r t ,x 1 ,x 2 )= f(r t ,x 1 ,x 2 ) + Af(x 1 ,x 2 ) (16) 

C ss( X l> X 2) = C ss( X l> X 2) + AC ss( X l> X 2)> ( 17 ) 

where f( r t> x i> x 2) and C ss (x 1 ,x 2 ) denote nomina i mo dels 

while A f( x v x i) an d AC ss (*i>*2) are the uncertainties, 
respectively. 

In the following, we will discuss those two uncertainties. 
Details are given below. 

1) Uncertainty in Compressor Torque: 

Based on the conservation of angular momentum in spool 
dynamics, the compressor load torque is required to overcome 
the duct wall friction [16]. For a compression system, Euler's 
turbine equation is used to describe the compressor torque 
which depends on the angular momentum of the fluid [18]. The 
schematic compressor load torque is adopted from [18] as 
shown in Fig. 5. As a result of turbo-machine driving at an 
angular speed, the compressor torque can be obtained by 
Euler's turbine equation as given by (e.g., [9], [18]): 



: <f>(r 2 v 02 - r^) , 



(18) 



where x and 2 are radii at the inducer tip and the inducer hub 

casing, respectively. In addition, the notation 0l denotes 
instantaneous tangential velocity of which the swirling fluid 

passes through the inducer tip at radius 1 while e2 is the 

tangential velocity at inducer hub casing with radius 2 . 
However, the various power losses might exist in the practical 
system such as incidence losses, fluid friction losses and inlet 
casing losses [9]. Those energy losses might procure the 
fluctuation behaviour into the compressor torque. Hence, such 
kinds of uncertainty and disturbance need to be considered in 
the dynamical model of a centrifugal compressor. 

In order to compensate the uncertainty in compressor 
torque, we introduce the sector nonlinearity method to guard 
the operating boundary of compressor torque. It is observed 
from Eq. (18) and Fig. 5 that the compressor's driving torque is 
associated with mass flow rate and spool angular speed. Based 
on the above discussion, we make the following assumption. 




Inducer tip 



Inducer hub casing 



Fig. 5 Compressor Torque inTurbo-Machine 

Assumption 1: Suppose that the compressor's driving 
torque with uncertainty satisfies the following condition: 



f(T t ,x 1 ,x 2 ) + Af(x 1 ,x 2 )\<M 1 \x 1 \+M 2 



(19) 



for all 1 and 2 in 
containing the origin. 



D e s Jt ^ where D denotes an open set 



M 1 \x 1 



M, 



Xj 



Note that, J "il A H and 
allowed local sector bounds for the 
torque function. 



are apparently the 
compressor's driving 



2) Uncertainly in Compressor Characteristic Maps: 

It is known that the stable operation region of the 
compressor characteristic can be collected from practical 
experiment. However, the unstable behaviour of compressor 
dynamics is hardly measured. Based on the structure of the 
centrifugal compressor, inexact unstable branch might cause 
the existence of uncertainty in the compressor characteristic 
maps. Thus, it is inevitable that such an uncertainty requires to 
be considered in the design of robust control scheme. For the 
centrifugal compression system, the characteristic map depends 
on the mass flow rate and spool angular speed. We have the 
next assumption for characteristic map. 

Assumptions 2: Suppose that the compressor characteristic 
map with uncertainty satisfies the following condition: 

| C ss (x 1 ,x 2 ) + AC ss (x 1; x 2 ) |< N, \ Xl | +JV 2 | x 2 | (20) 

for all x x and x 2 in D e $R 2 , where D denotes an open set 
containing the origin. 



Note that, N 1 \x 1 \ and JV 2 



in Assumption 1 above 



denote the local sector bounds for compressor characteristic 
map. 

B. Robust Control Design 

To fulfill the stabilization design for the centrifugal 
compressor dynamics, we employ the robust control scheme to 
carry out the suppression of surge phenomenon. According to 
Assumptions 1 and 2, the nonlinear terms of the system with 
the appearance of uncertainties in both of compressor's driving 
torque and characteristic map are assumed to be limited in local 
sector bounds. Furthermore, it is clear that the system (7)-(9) is 
controllable. In this paper, we employ the local sector approach 
to design the robust controller for the system (7)-(9) with both 

driving torque 1 and close-couple valve 2 . 



Let 



V(X): 



2"a ' 



(— *r + 



1 2 1 2x 
= *2+-*3) 

b c 



(21) 



be a Lyapunov function candidate for system (7)-(9). The time 
derivative of the Lyapunov function candidate is given by 



1 1 

V(x) = —x 1 x l 
a 



1 



\ 2 A 2 



^3 A 3 



b - c 
-- xj f(r te ,x lt x 2 ) + Af{x 1 + x 2 ) + uj 
+ x 2l c ss( x i> x i) + AC ss (x 1; x 2 ) + u 2 ] 

+ x -i\F(re>Ve)- F (re> x 3+¥e)\ 

< MjX 2 + x l u l + (M 2 + JV-l) I x 1 1 • | x 

+x 2 u 2 + xJ_F(y e ,y/ e ) - F{y e ,x 3 + y/ e )] 



2 I +JV 2*2 



(22) 
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Assume F is a strictly increasing function with respect to 3 . 
It is obvious to have that x d F (r e .¥e) - F(y e> x 3 + Ve )} < for 

ail x 3 whUe x 3 [F(r e , ¥e )-F(r e ,x 3 + ¥e )-] = o hM at x 3 = 

[12]. Let the control inputs 1 and 2 be chosen as 



The notations " and '< , != l»---» 5 are non-dimensional 
parameters, and a denotes the slip factor. In the modeling of 
the centrifugal compressor, the characteristic map given in Eq. 

(29) holds at " . For the case of deep surge involving 
negative mass flow, the compressor characteristic of deep 



«i = ~Vi 



and 



Li n — Kn At 



(23) 
(24) 



surge for ■ 
where 



■<0 



is shown as ([7]-[8]) 



For 1 and 2 . We then have 

V(x) < (M 1 - k,)x 2 + (M 2 + N 1 )\x 1 \-\x 2 \ +(N 2 - k 2 )x 2 2 



y/ c0 {a>,f) = (l+arj c0 do} 2 y c - 1 -1 



(30) 



(31) 



and 



= (M 1 -k 1 )[|x 1 | + 



(M 2 + JVQ- 1 x. 
2(M 1 -fc 1 ) 



2 h2 



2a 



+ [N ? -k 



(m 2 + N,y 



1 A.T JAl 

- - 4(Mi-K) 



(25) 



V(x) 



It is clear fromthe above discussions that *• ' in (25) will 

k k 

be a negative definite function, while both 1 and 2 are large 

enough. We then have the next results. 

Theorem 1: Suppose that Assumptions 1 and 2 hold. The 

equilibrium points Xe of system (7)-(9) will be locally 

asymptotically stabilized by a linear control law given in (23) 

k k 

and (24) with l and 2 being large enough. 

Note that, the proposed control scheme given in Theorem 1 

can be achieved without having exact knowledge of the 

C (x x ) 
compressor characteristic map ss ^ v 2J and compressor 

torque function ' ^ Tt ' x ' 2 > . However, it does require to have 
the sector bound to cover the local behaviour of nonlinear 
terms with the uncertainties. 

IV. Numerical Simulations 
In order to check the effectiveness of the robust controller 
design on system performance, we adopt the centrifugal 
compressor model with same parameter values from [7]-[8] for 
numerical study. Here, the compressor torque is given by 



r(&>,^) = a\0\co 
and throttle function 

F{y,y/) = yjyr . 

In addition, the compressor characteristic maps ^ c ^ ' 
given by 

y/ c (a>, <t>) = [1 + at] c (a>, (/>)da> 2 ] ^ - 1 , 
where 



(26) 



(27) 



(28) 



ti c (m,</>)-- 



GO) 



am 2 + \(f l( o- f 2 <f>) 2 +U<jw- U) 2 +<f> 2 (h + f 5 ) ' 



(29) 



a +2a+ f t 
The value of parameters are given as: 

(a, b, c, d) = (6.97, 310.81, 9.37, 0.38) 
(/i. U, U> U, U) = (0.44,1.07, 2.18, 0.17, 0.12) 
7c =1.4, p = 3 and a = 0.9 

To compare with the results shown in Fig. 3, for the 
following numerical simulation we choose ^~ ' with 

• w i a;- (5 n ,^ n ,^n) = (0.5, 0.345, 0.22) 
initial condition v O'VO'ru^ <. > > i 

A. Open-Loop Dynamics 

Based on the Fig. 3, we choose a different equilibrium point 

to simulate the dynamical behaviour of the centrifugal 

compressor. The numerical simulations of timing responses 

and phase plane for the open-loop system are obtained by using 

code Matlab and are shown in Figs. 6-8. As depicted in Fig. 6, 

t , „ .,.- • (<sL«LiO = (0.766608, 

the system equilibrium *• <""' rE ' *• 

0.289736, 0.846) , t. =0.1999. t a* k , *• n 

y for f is found to be asymptotically 

stable. The surge instabilities are observed in Figs. 7 and 8. In 

Fig. 7, we do observe that system behaviour reveals period-two 

oscillation for the system equilibrium 

(£ e ,&,^ e ) = (0.606924,0.21964,0.4864) when t t = 0.12 ^ 

addition, the multi-period oscillation of centrifugal compressor 
dynamics is obtained for the equilibrium 

point (S e ,&,V e ) = (0.665538, 0.244524, 0.6026) when 

f ~ ' . The corresponding timing responses and phase 

portrait are shown in Fig. 8. 

Under above analysis of the uncontrolled centrifugal 
compressor model, the system might fall into a violent surge 
phenomenon with respect to the variation of driving torque 
and/or throttle opening. It can be observed from Figs. 6-8 that 
slight variation of the driving torque and/or the throttle opening 
might cause the surge dynamics to occur in the centrifugal 
compressor system. 

B. Closed-Loop Dynamics 

In the following, the robust control schemes developed in 
Section III are employed to stabilize the system behaviour 
which considers the system uncertainty in the centrifugal 
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compressor dynamics. The effectiveness of the proposed robust 
control law will also be demonstrated. 

Consider the system (4)-(6) with control inputs > and 2 
as given by 

co = a[T [ -T{a),<P) + u 1 ] (32) 

<j> = b[if/ c (co,0)- l t / + U 2 ] (33) 

W = c[</>-F(y, V )]. (34) 

Here, the control inputs 1 and 2 denote driving torque 
and close-couple valve, respectively. According to Theorem 1, 
we have the control law as given by 



compressor's driving torque and characteristic map. The 
system (7)-(9) can then be rewritten as 



u^ =-k l {co-co e ) 



and 



u 2 = -k 2 {<f> - e ) , 



(35) 



(36) 



«K 



where e and Ve denote the value of the corresponding 
unstable equilibrium point at which system might exhibit 
oscillative behaviour. 

1) Controller Design without System Uncertainty: 

First, we adopt the centrifugal compressor model without 

considering system uncertainty. The control inputs ! and 

2 given in (35) and (36) are used for system(32)-(34) without 
considering system uncertainty. It is clear to have the bounds 

for the functions f^t^i^i) ana c ss ( x i> x 2) t0 meet 

Assumptions 1 and 2. For the control design, we choose 

fc = 0.44 A kj = 0.59 T, t t , . , 

1 and 2 . The parameters for the numerical 

study are obtained for three different equilibrium 

points (® e .^.^) = (0-557403, 0.199432, 0.4008) _ (5 e .£ f 

y/ e ) = (0.606924, 0.219694, 0.4864) and (& e ,0 e ,y/ e ) = 

(0.665538,0.244524,0.6026) . , . t , . c . Q „ ., 
v ' ' ' . As depicted in Figs. 9-11, the 

timing trajectories of system states are pushed toward the 

system equilibrium by the control inputs 1 and 2 . As 
depicted in Figs. 4 and 9, the system state is clearly pushed 
from period-one oscillation to a stable operating point by 

control inputs l and 2 . Compared with the surge dynamics 
depicted in Fig. 7, it can be observed that the period-two 
behaviour is stabilized to the equilibrium point as shown in Fig. 
10. Similarly, the complicated dynamics depicted in Fig. 8 is 
also stabilized by the combination of the driving torque and 
close-coupled valve control as shown in Fig. 11. In addition, as 
depicted in Fig. 12, we consider that the controller is switched 

at time f = ^ . It is observed from Fig. 12 that the system 
instabilities prior to the control action are diminished after 
adding control effort. 

2) Controller Design with System Uncertainty: 

Next, we present the numerical results of the system (7)-(9) 
with the robust control law. Following the control design 

procedure proposed in Section III, the control inputs 1 and 2 
are used to compensate the system uncertainties embedded in 



x 1 = a[f(T t ,x 1 ,x 2 ) + Af(x 1 ,x 2 ) + u l ] 

x 2 = b[C ss (x v x 2 ) + AC SS (x v x 2 )-x 3 +u 2 ] 

x 3 =c[x 2 +g{y e ,x 3 )] 



(37) 
(38) 
(39) 

Motivated by [12], the system uncertainties are chosen as 
Af(x 1 ,x 2 ) = 0.1sin(10x 2 ) (40) 

AC ss (x 1 ,x 2 ) = 0.1sin(10x 2 ). (41) 

The characteristic maps of centrifugal compressor with 
system uncertainty are shown in Fig. 13. In those figures, the 
dot line denotes the characteristic maps with uncertainty and 
the solid line is the nominal one. It is not difficult to find the 

bounds for the functions f ( T r> x i' x 2) + A f ( x i> x i) ana 
C ss (x 1 ,x 2 ) + AC ss (x 1) x 2 ) tQ met the Assumptions x and 2 . 



First, we use the control gain 



Jq = 0.44 



and k 2 



0.59 



to 



stabilize the spool dynamics of the system (37)-(39) with the 
addition of system uncertainties. Numerical simulations are 
shown in Fig. 14. It is observed from Fig. 14 that those control 
gains can not compensate the system uncertainties for the 

,., . . , (w e ,S.,w.) = (0.665538, 0.244524, 0.6026) 

equilibrium point v e>re>re/ v > > > 

The system dynamics still exhibits oscillative behaviour. Based 

on the bounds for the uncertainties (40) and (41), we choose 

k = 5 
the control gain to be larger. That is, to let x ' and 

k =13 
2 ' . The simulation results depited in Fig. 15 show that 

the oscillative dynamics is stabilized by using the larger gains 

k k 

of control 1 and 2 . As shown in Fig. 16, a switching control 

scheme is added to demonstrate the effectiveness of the 

proposed robust control law. The control will be applied to 

system at time f = ^ . It is observed that the system instabilities 
can be diminished by robust controller. The numerical 
simulations are found to agree with the analytical results of 
Theorem 1. 
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Fig. 6 Dynamical Behaviour in Open-Loop System for Equilibrium Point 
(a> e ,&, iO = (0.766608, 0.289736, 0.846) 
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Fig. 7 Dynamical Behaviour in Open-Loop SystemfFor Equilibrium Point 
(w e ,^,yf e ) = (0.606924, 0.2197, 0.4864) 
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Fig. 8 Dynamical Behaviour in Open-Loop System forEquilibriumPoint 
H , $ c , y/ e ) = (0.665538, 0.244524, 0.6026) 
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Fig. 9 Dynamical Behaviour in Closed- Loop System for Equilibrium Point 
H , tj> e , W e ) = (0.557403, 0. 199432, 0.4008) 
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Fig. 10 Dynamical Behaviour in Closed-Loop System forEquilibriumPoint 
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Fig. 13 Compare withthe Uncertainty in the Characteristic Maps 
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Fig. 14 Closed-Loop System with Uncertainties for Equilibrium Point 
K ,fc,y/ t ) = (0.665538, 0.244524, 0.6026) 
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Fig. 15 Dynamical Behaviour with Robust Control Scheme for Equilibrium 
Point (a> e ,<j> e ,y/ e ) = (0.665538, 0.244524, 0.6026) 



Fig. 16 Dynamical Behaviour with Robust Control Scheme for Equilibrium 
Point (a e ,^„,y/ e ) = (0.665538, 0.244524, 0.6026) 

Conclusions 

A robust control design is proposed in this paper for the 
control of surge phenomenon in centrifugal compressor 
dynamics via local sector nonlinearity approach. The proposed 
control scheme was shown to be robust with the appearance of 
system uncertainties embedded in the unstable operation region 
of centrifugal compressor. In addition, the proposed robust 
controller design can achieve the system stabilization without 
having exact knowledge of characteristic map function and 
compressor's driving torque function. Based on such an 
advantage, the proposed control design might provide an 
implementable strategy in the practical applications. 
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